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The proposed semiclassical theory predicts two types of oscillations in the flow of current injected 
from a point source near a ballistic p-n junction in graphene in a strong magnetic field. One 
originates from the classical effect of bunching of cyclotron orbits of electrons passing back and 
forth across the p-n interface, which displays a pronounced dependence on the commensurability 
between the cyclotron radii in the n- and p-regions. The other effect is caused by the interference 
of monochromatic electron waves in p-n junctions with equal carrier densities on the two sides and 
it consists in magneto-oscillations in the current transmission through the interface with periodicity 
similar to Shubnikov-de Haas oscillations. 

PACS numbers: 72.80.Vp, 73.40.Lq, 73.43.Qt 




Graphene is a gapless semiconductor with charge car- 
riers that behave like massless Dirac particles [TJ [2] , in 
which it is possible to locally control the carrier density 
and type using local electrostatic gates, and create p-n 
junctions [3H6]. Due to the Dirac- type properties of the 
carriers, such p-n junctions are highly transparent for in- 
coming electrons [TJ [8] . Also, it has been suggested that 
p-n junctions in ballistic graphene may be able to focus 
flow of electrons injected through a point-like source [9], 
which is the result of the inverted dispersions e c = vp 
and e v = U — vp, for electrons in the conduction and 
valence bands on opposite sides of the junction, and con- 
sequently, an effectively negative index of refraction for 
electron trajectories. 

In this paper we study the features in the flow of elec- 
trons injected in the vicinity of a p-n junction in ballistic 
graphene in a strong magnetic field, resulting from the 
negative refraction of electrons crossing the p-n inter- 
face. Reflected and transmitted electrons follow trajec- 
tories which combine elements of skipping orbits [TQHT2] 
(formed by repeatedly reflected electrons) and snaking 
orbits (formed by electrons repeatedly crossing the p-n 
interface). In a recent experimental work [13] , signatures 
of such snaking orbits were observed in transport mea- 
surements on gated graphene p-n junctions. Also, It has 
been noticed [TIJ [15] that when skipping orbits originate 
from a point-like source, they bunch into caustics and 
exhibit focusing of charge flow at periodically repeated 
cusps. Here, we show that periodic appearance of caus- 
tics and their cusps is also characteristic for the combined 
skipping-snaking orbits of electrons propagating along a 
p-n interface in ballistic graphene, and we analyse how 
the resulting singularities in the flow of current injected 
using a point-like source depend on the relative value of 
the densities of carriers (electrons/holes) on the two sides 
of the p-n junction. A typical distribution of current cal- 
culated in such a system is illustrated in Fig. [I] where 
the formation of the p-n junction is defined using a pair 
of split gates providing a potential step U = v(p n + p p ) 




FIG. 1. Current distribution in graphene p-n junction in a 
perpendicular magnetic field, with electrons injected isotrop- 
ically at the Fermi energy from the point source S. 



for electrons in graphene, and the maxima in the current 
intensity are marked as the bright spots. Also, we find 
that for a junction with equal carrier densities on the 
two sides, the distribution of current displays additional 
quantum magneto-oscillations. 

The proposed theory is based on the analysis of the 
families of orbits of electrons injected, e.g., on the n- 
side of the p-n junction at a distance xq from it. The 
radii of the electron cyclotron orbits on either side of the 
interface are given by r n ( p ) = p n {p)/eB, where p n (p) are 
the electron Fermi momenta in the doped regions. Also, 
as in Ref. [8 , we assume the p-n junction potential step U 
to be sharp on the scale of the electron Fermi wavelength. 
For each individual trajectory, such as shown in Fig. [2j 
an electron leaves the source at a certain angle < 
$ < 2ir. The following path of the electron consists of 
a sequence of circular segments on either the n-doped 
(with radius r n ) or the p-doped (with r p ) side, which are 
matched at the interface by the kinematically prescribed 
Snell's law with negative refraction index [9], 
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Here, the values of electron momenta p n ( p ) (where 
v(p n + Pp) = U) are absorbed in the values of the cy- 
clotron radii r n ( p ), and is the angle of incidence and 
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FIG. 2. Branching of trajectories of electrons propagating in 
a magnetic field along an asymmetric, r n ^ r p (top), and 
symmetric, r n = r p (bottom) p-n junction. 



specular reflection, while 0' is the angle of refraction. The 
sign of the refraction index in Eq.Q is negative because 
the group velocity is v c = vp/p in the conduction band 
and v v = — vp/p in the valence band. The latter relation 
also prescribes the reversal of the direction of the elec- 



tron's angular velocity as it crosses the interface, making 
electrons drift along the p-n interface in a magnetic field 
(Fig.§. 

To describe the distribution of the current carried by 
the drifting electrons, one has to take into account that 
electron trajectories branch at the interface, with the 
probability Wt(0) to be transmitted and Wr = 1 — Wt 
to be reflected. For electrons arriving at a sharp p-n 
interface with angle of incidence |0| < C [3|H], 
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and Wt = for \0\ > (9 C , where smO c = 1/n and k = 
r n /r p . Branching of the electron trajectories can be then 
labeled by a sequence, 



M = {si, s 2 , , sjv(A4)-i? s n(m)}i 
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where si = +1 if the i th circle segment is on the n-side 
and Si = —1 if it is on the p-side of the junction. Together 
with the sequence uniquely defines the path with N 
circle segments, as shown in Fig. [5] Using Snell's law 
and elementary geometry, one can show that the centers 
of all cyclotron orbits whose segments belong to the se- 
quence hA are positioned at the same distance X{d) from 
the p-n interface. 

To calculate the density g(x, y) of trajectories on the 
graphene sheet, we define a generating function [15], 



F(x, y, ti, M) = [x- X(V)} 2 + [y- Y(0, M)\ - r 2 (M), 
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X(ti) = x + r n smti, r(M)=s N(M) r -^j^ + r -^^, Zi (M) = Si (M) + a i+1 (M). (4) 
For each D and M, a point (x, y) lies on the N th segment of the path if F = 0, thus, the density of trajectories is 
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where we take into account all sequences whose final seg- 
ments pass through (x, y) for all possible which solve 
equation F = 0, and tir^t) is the number of reflec- 
tions (transmissions) contained in the sequence A4. 

Figure [3] shows the spatial distribution of q(x, y), calcu- 
lated from Eq. (|6|, which illustrates spreading of current 



flow injected from a point-like source along the p-n inter- 
face, for various combintations of system parameters: Xq, 
r n and r p . It reflects the periodic appearance of singular- 
ities in the skipping-snaking orbits (caustics and caustic 
cusps), weighted with the reflection/transmission proba- 
bilities. In all panels in Fig.|3j caustics are seen as bright 






FIG. 3. Calculated density of current injected from a point-like source (bright spot) near the p-n junction. Electrons are 
injected isotropically at a distance xo from the interface and travel along the interface undergoing transmissions and reflections. 
Spatial oscillations of current distribution are caused by a periodic appearance of caustics and cusp singularities 



lines, indicating the points (x c , y c ) on the 2D plane where 
both F = and dF/dd = 0: 
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a(&) = dX{ti)/dd, P(ti,M) = dY(&,M)/dti. 

Cusps, representing the bright spots of "magnetic focus- 
ing" of the electron flow [TJJ [15], appear at the points 
where, additionally, d 2 F/d 2 d = 0. 

When r p = (1 st column of Fig. [3]), all electrons strik- 
ing the interface on the n-side are reflected, so that, 
classically, current would flow only on the side of the 
source, displaying singularities (single cusps as well as 
cusp triplets) in the distribution specific for the families 
of simple skipping orbits [15]. As r p is increased (2 nd col- 
umn), electrons start getting transmitted to the p-side, 
leading to caustics and cusp-like features on both sides 
of the interface. As r p is made comparable to r n , the 
critical angle 6 C disappears, and the strength of features 
produced by reflected electrons diminishes significantly 
(3 rd column). In fact, for xo = r n = r p , the caustics and 
cusps produced by skipping trajectories vanish, and we 
only find the features produced by snaking orbits, lead- 
ing to a doubling of the period of appearance of cusps, 



as compared to the case of purely skipping trajectories. 
As r p is increased further (4 th column), the period of 
appearance of cusps of snaking orbits and the number 
of cusps of skipping orbits between consecutive cusps of 
snaking orbits increase. Note that, in the general case 
of r n ^ r p , the trajectories for different paths are not 
commensurate for all angles. This determines a complex 
pattern of current distribution. In special cases where 
the radii r n and r p are commensurate (r p = mr n ), the 
singularities in the current distribution on the side of the 
source or the opposite side of the p-n junction display 
long range periodicity, at mr n . The evolution of the cur- 
rent distribution upon the change of magnetic field and 
carrier density is shown in Fig. |4j for current measured 
at a distance of 4x from the source. 

In the special case of a symmetric p-n junction, with 
r n = r p , it is also possible to develop a coherent semi- 
classical description of current distribution, taking into 
account quantum interference effects. This involves tak- 
ing into account the phases acquired by electrons upon 
transmission/reflection at the interface and during free 
propagation of Dirac electrons. As shown in Fig. |2j for 
the case r n = r p , all families M of electron trajecto- 
ries for any xq and d belong to a periodic sequence of 
circles, which enables one to describe the wave propagat- 
ing along such paths using a pair of amplitudes, ip n (j^) 
and ip p (j,$) for an electron on the n- or p-side of the 



4 




FIG. 4. Calculated local density of current at the point 
(xo,4xo), located directly below the source at a distance of 
4xo from it. The magnetic field B and carrier density n n ( p ) 
dependence of the local current density are encoded in the 
dependence of r p /r n = \fn p /n n and r n /xo oc 1/B. 



interface after j — 1 encounters with it. The wave am- 
plitudes i/j n ( p )(j, #) on the j th circle segment on the n(p) 
side, immediately after the (j — l) th interaction with the 
interface, were found as, 
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Here the scattering matrix S describes the transmis- 
sion and reflection of Dirac spinors at the interface, and 
S their propagation between interactions with the in- 
terface. These matrices were calculated for the plane 
wave Dirac spinors of electrons in graphene [8] , and 
hence the product SS includes the electron Berry phase. 
The phase </>o is acquired by electrons while propagat- 
ing between the source and the interface, and n ( p ) = 

J \—eA/ h + p n (p) I h) ' dl are phases acquired by electrons 
propagating in a magnetic field along a single circular 
segment. As segments on the p- and n-sides always form 
full circles, Fig. [2J we find that <p n — <p p = 7rp 2 /(HeB). 

Due to the above-discussed feature, the drift of the 
Dirac electron along the p-n interface described by 
Eq. ([7]) has the following interesting property: for any 
$ and #o, both the matrix V = SSE and the matrix 
W = [SE] 2 are diagonal. According to Eq. (H), af- 



ter 2N encounters with the interface, the amplitudes 
ijj n ( p }(2N + of the Dirac electron on the n(p)-sides 
are determined by the evolution matrix (S6') 2Ar_1 S = 
VW n-i Therefore, if p 2 /(heB) = (21 + 1), the elec- 
trons will appear only on the same side of the junction 
as their source after any even number of encounters with 
the p-n interface. As a result, we predict rapid quantum 
oscillations of the electron current between the two sides 
of the junction. These oscillations develop on a scale 
linear in the inverse magnetic field with the same 

period A(B~ 1 ) = 2he/p 2 as Shubnikov-de Haas oscilla- 
tions [2]: the phase difference, (j) n — (j) p , coincides exactly 
with the phase acquired by an electron with momentum 
p after it has propagated along a closed cyclotron trajec- 
tory in the bulk of the graphene flake with carrier density 
n = p 2 /(7rh 2 ). 

At a finite temperature T > fw;/(4/c#r), r = r n = r p , 
thermal smearing of the Fermi level results in the loss of 
monochromaticity of the electron source, which attenu- 
ates the above-mentioned intereference effect, leading to 
the classical results applicable to the incoherent propaga- 
tion of electrons. Then, incoherent propogation gives rise 
to a classical density profile for the distribution of elec- 
trons injected by the DC source, as defined in Eq. §6§. For 
the case of r n = r p , summation of the series in Eq. ^ 
can be performed using a recursive method similar to 
Eq. ([7]), resulting in 
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where $ is chosen to solve F = for the point (x,y), 
and which produces an alternating pattern of weak and 
strong cusps, as shown in Fig. [3] 

In conclusion, we predict periodic spatial modulation 
of current injected from a point source near a p-n junc- 
tion in graphene in a magnetic field, which originates 
from caustic bunching of skipping-snaking orbits. Ex- 
perimentally, such oscillations would appear as magneto- 
oscillations of conductance of a two-terminal device with 
two point contacts. We also find that, for the case of com- 
mensurate size of electron cyclotron orbits in the n-type 
and p-type regions, the current distribution undergoes 
rapid quantum oscillations between the two sides. 
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